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Abstract. We present the performances of a 330 g zinc molybdate (ZnMoOzt) crystal working as scin- 
tillating bolometer as a possible candidate for a next generation experiment to search for neutrinoless 
double beta decay of 100 Mo. The energy resolution, evaluated at the 2615 keV 7-line of 208 T1, is 6.3 keV 
FWHM. The internal radioactive contaminations of the ZnMo04 were evaluated as < 6 /iBq/kg ( 228 Th) 
and 27±6 ^iBq/kg ( 226 Ra). We also present the results of the a vs /3/7 discrimination, obtained through 
the scintillation light as well as through the study of the shape of the thermal signal alone. 

PACS. 23.40.Bw Weak interactions in /3 decay - 29.40.Mc Scintillation detectors - 07.57.Kp Bolometers 



1 Introduction 

The neutrinoless Double Beta Decay (0VDBD) is a nuclear 
process that, if observed, would establish that the total 
lepton number is not a conserved quantity and that the 
neutrino is a Majorana particle, and would set the abso- 
lute mass scale of neutrinos [T|[2][3] • Plenty of experiments 
are now in the construction phase and many others are in 
R&D phase g]. Very recently EXO [5] and KamLAND- 
Zen [B] set very competitive limits on the Oz/DBD half-life 
of 136 Xe. 

The main challenges for all the different experimental 
techniques are the same [7]: i) increase of the active mass, 
ii) decrease of the background, and iii) improvement of 
the energy resolution. 

Thermal bolometers are ideal detectors for this kind of 
research: crystals can be grown with a variety of interest- 
ing DBD-emitters and multi-kg detectors can be operated 
with excellent energy resolution [B] which, perhaps, repre- 
sents one of the most critical aspects for next generation 
experiments. 
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The Cuoricino experiment [3] searched for the 0VDBD 
of 130 Te operating 62 Te02 bolometers. The Cuoricino 
data demonstrated that the background in the region of 
interest is dominated by radioactive contaminations on 
the surfaces facing the detectors, a particles produced by 
these contaminants can lose a fraction of their energy in 
the host material, and the rest in the detector, thus pro- 
ducing a flat background from the energy of the decay 
(several MeV) down to the 0VDBD region [TO]. More- 
over simulations show that this contribution will largely 
dominate the expected background of the CUORE exper- 
iment [TTl[T2"] in the region of interest, since Te02 bolome- 
ters do not allow to distinguish a particles from the elec- 
trons emitted in the fVDBD. 

The natural way to discriminate this background, is 
to use scintillating bolometers [13] ■ In such devices the 
simultaneous and independent readout of the heat and the 
scintillation light signals permits to discriminate events 
due to /3/7, a and neutrons interactions thanks to their 
different scintillation yield. 

100 Mo is a very interesting /3/3-isotope because of its 
large transition energy Qpp = 3034 keV and a consider- 
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able natural isotopic abundance S = 9.67%. Several in- 
organic scintillators containing molybdenum were devel- 
oped in the last years. Among them, ZnMo04 was re- 
cently grown [T3] and the first cryogenic detector gave 
very promising results [T5] . Some scintillating crystals, the 
molybdates in particular, show a very peculiar feature: the 
thermal pulse induced by an a particle shows a slightly 
faster decay time with respect to the one induced by 7 
interactions [TB]. This feature seems to be explained [IT] 
by the relatively long scintillation decay time (of the or- 
der of hundreds of fis) observed in some scintillating crys- 
tals. This long decay, combined with a high percentage 
of non-radiative de-excitation of the scintillation channel, 
will transfer phonons (i.e. heat) to the crystal. This ex- 
tremely tiny, but measurable, time dependent phonon re- 
lease has a different absolute value for isoenergetic a and 
/3/7 particles due to their different scintillation yield. 

It was very recently measured that the a vs 7 sepa- 
ration on a 29 g ZnMo04 crystal can reach an extremely 
high efficiency |19j using the pulse shape discrimination 
(PSD) alone, while the separation based on the scintilla- 
tion light, even on smaller crystals (5 g), shows a smaller 
efficiency [2"0] . 

A next generation experiment, nonetheless, will need 
to run detectors with a considerably larger mass (of the 
order of few hundreds of grams each) . It is not straightfor- 
ward to foresee the performance of the particle discrimi- 
nation method over a significantly larger sample. This is 
due to the fact that the PSD is sensitive to the signal to 
noise ratio, and a larger mass of the absorber crystal leads 
to a smaller signal amplitude (oc (detector mass) -1 ). The 
same holds (in principle) for the discrimination based on 
the scintillation light: usually (especially in the case of non 
transparent crystals) the larger the crystal, the smaller the 
light output. Moreover one has to consider that this com- 
pound is characterized by an extremely tiny Light Yield 
(LY): 1+2 keV/MeV. 

The purpose of this work is to study the most im- 
portant parameters (energy resolution, a vs 7 discrimina- 
tion, internal radiopurity) on a crystal whose size (330 g) 
matches very closely the requirement of a next generation 
O^DBD experiment. 

2 Experimental set-up 

The 330 g ZnMoC>4 crystal studied in this work was grown 
in the Nikolaev Institute of Inorganic Chemistry (NIIC, 
Novosibirsk, Russia). Starting material for the crystal gro- 
wth were high purity ZnO (produced by Umicore) and 
M0O3, synthesized by NIIC. 

Crystals are grown by the low-thermal-gradient Czo- 
chralski technique (LTG Cz) [TH] from a melt contained in 
a 80 mm diameter platinum crucible. The sample was cut 
from one of the first large-size crystals grown by the LTG 
Cz. The development of the growth process of ZnMo04 
is at the very beginning, and the shape of the crystal is 
non-optimal. The final form of the crystal was chosen as 
a compromise between a large crystal size and the mini- 
mization of visible defects. The shape of the crystal sample 




Fig. 1. Set-up of the detectors. The ball-bonded Au wires 
are crimped into "male" Cu tubes (pins) and inserted into 
ground-insulated "female" Cu tubes. Custom wires from de- 
tectors towards cryostat are not drawn. A section of the light 
detector and of the reflecting sheet is not drawn for a better 
understanding. 



used in the current work is an irregular polyhedron with 
5 sides whose cross section can be roughly assimilated to 
a 45x45 mm 2 square. The only parallel faces are the top 
and bottom ones. All the surfaces are matted, except the 
one faced to the light detector, that was polished. Unlike 
previous small samples (T9"l[2"0] that were colorless, this 
crystal shows an uniform orange tint. 

The ZnMo04 crystal is held by means of four S-shaped 
PTFE supports fixed to cylindrical Cu columns: two on 
the top and two on the bottom. The crystal is surrounded 
laterally and on the bottom part (with no direct thermal 
contact) by a plastic reflecting sheet (3M VM2002). The 
temperature sensors are 3x3x1 mm 3 Neutron Transmuta- 
tion Doped (NTD) germanium thermistors, the same used 
in the Cuoricino experiment. For redundancy we decided 
to use two thermistors. Each of them is thermally coupled 
to the crystal via 9 glue spots of « 0.6 mm diameter and 
ss 50 /im height. 

At the working temperature of our bolometers (10-^30 
mK), no "standard" light detectors can work properly. 
The best way to overcome this problem is to use a second 
-very sensitive- "dark" bolometer that absorbs the scin- 
tillation light giving rise to a measurable increase of its 
temperature [21] . Our Light Detector (LD) consists of a 
50 mm diameter, 260 /jm thick pure Ge crystal absorber 
facing the polished surface of the crystal. A schematic view 
of our set-up is presented in Fig. [TJ 

The detectors were operated deep underground in the 
Gran Sasso National Laboratories in the CUORE R&D 
test cryostat. The details of the electronics and the cryo- 
genic facility can be found elsewhere [2"2]|23(I24"]. 
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The heat and light pulses, produced by a particle inter- 
acting in the absorber and transduced in a voltage pulse 
by the NTD thermistors, are amplified and fed into an 
18-bits NI-6284 PXI ADC unit. The trigger is software 
generated on each thermistor and when it fires 1 s long 
waveforms, sampled at 2 kHz, are then saved on disk. The 
time window for the LD is shorter and corresponds to 250 
ms. Moreover, when the trigger of a ZnMo04 thermistor 
fires, the corresponding waveform from the LD is recorded, 
irrespective of its trigger. 

As one of the main goals of the measurements was to 
test the a vs discrimination capability of this large 
crystal, a 238 u/ 234 U a source was faced to the crystals, on 
the opposite side with respect to the LD. The source was 
covered with a 12 /xm thick polyethylene film, in order to 
smear the a's energies down to the 100 Mo Q^-value. Since 
mounted close to the detector, this source is responsible 
for an increase of the a background that could in principle 
spoil our sensitivity to the intrinsic contamination of the 
ZnMoC>4 crystal. However, the straggling of alpha parti- 
cles inside the polyethylene film "shifts" the 2 38tj/234tj 
a-particles toward lower energies (i.e. below 4 MeV) re- 
moving a possible interference with the 4 MeV peak that 
should appear in the case of a 232 Th bulk contamination 
of the crystal. 

The 7 calibration of the ZnMoC>4 crystal is performed 
through removable 228 Th and 40 K sources inserted be- 
tween the dewar housing the cryostat and the external 
lead shield. The energy calibration of the LD is achieved 
thanks to a permanent 55 Fe X-ray source, producing two 
X-rays at 5.9 and 6.5 keV, faced closely to the LD. 

2.1 Data Analysis 

The amplitude and the shape of the voltage pulse is de- 
termined by the off-line analysis that makes use of the 
Optimum Filter technique [2"51l2l)] . The signal amplitudes 
are computed as the maximum of the filtered pulse. The 
amplitude of the light signal is estimated from the value of 
the filtered waveform at a fixed time delay with respect to 
the signal of the ZnMoC>4 bolometer, as described in de- 
tail in Ref. [37] . The signal shape is evaluated on the basis 
of four different parameters: tr, Td, TVL and TVR. tr 
(the rise time) and tr> (the decay time) are evaluated on 
the raw pulse as (tgo%-tio%) and (t3o%-tgo%) respectively. 
TVR (Test Value Right) and TVL (Test Value Left) are 
computed on the filtered pulse as the least square differ- 
ences with respect to the filtered response function^ of the 
detector: TVR on the right and TVL on the left side of 
the optimally filtered pulse maximum. These two param- 
eters do not have a direct physical meaning, however they 
are extremely sensitive (even in noisy conditions) to any 
difference between the shape of the analyzed pulse and 
the response function. The detector performances are re- 

1 The response function of the detector, i.e. the shape of 
a pulse in absence of noise, is estimated from the average of 
a large number of raw pulses. It is also used, together with 
the measured noise power spectrum, to construct the transfer 
function of the Optimum Filter. 



Table 1. Technical details for the ZnMoCU bolometer (Ther- 
mistor 1 and Thermistor 2) and for the LD. Signal represents 
the absolute voltage drop across the thermistor for a unitary 
energy deposition. 



Crystal 


Signal 


FWHM base 


TR 


TD 




{/iV/MeV] 


[keV] 


[ms] 


[ms] 


ZnMo0 4 -l 


16 


3.6 


12.8 


59.8 


ZnMo0 4 -2 


14 


3.7 


12.0 


60.3 


LD 


1800 


0.20 


3.2 


8.2 




Energy [keV] 

Fig. 2. Calibration spectrum obtained by exposing the 
ZnMo0 4 crystal to the 228 Th source for 80 h. The peak at 
2615 keV of 208 T1, magnified in the inset, shows a FWHM res- 
olution of 6.3 keV. 



ported in Tab. [TJ The baseline resolution, FWHMj, nse , is 
governed by the noise fluctuation at the filter output, and 
does not depend on the absolute pulse amplitude. 

As mentioned above, the use of two thermistors on 
the same absorber is often made for redundancy. In this 
case we took advantage of the similar performance of both 
of them, using their sum. This technique is useful in the 
case that the noise fluctuations of the two thermistors are 
not correlated, meaning that these fluctuations are not ac- 
tual temperature fluctuation of the crystal. This technique 
can be used in two different ways. One can off-line com- 
bine the energies measured by the two thermistors into 
a "weighted" energy estimator [28] , linear combination of 
the two thermistors. Or, as in this case, one can sum the 
two signals at hardware level and treat the obtained signal 
as an additional independent channel. In our case the sum 
is performed after the two signals are amplified and just 
before they are fed into the acquisition. The calibration 
spectrum obtained on the sum of the two thermistors is 
presented in Fig. [5] 

The baseline energy resolution, FWHMf, ase evaluated 
on the sum (ZnMo04-Sum) is 2.6 keV, slightly better with 
respect to the ones reported in Tab. [I] The FWHM energy 
resolutions obtained at different energies are reported in 
Tab.H 
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Table 2. FWHM energy resolutions of the ZnMo0 4 detecto: 
evaluated on the two thermistors and on their sum. 





ZnMo0 4 -l 


ZnMo0 4 -2 


ZnMo0 4 -Sum 




[keV] 


[keV] 


[keV] 


583 keV 


4.1±0.7 


3.0±0.5 


2.9±0.4 


911 keV 


4.9±0.4 


4.7±0.5 


4.0±0.4 


1461 keV 


4.9±1.5 


5.4±1.2 


4.9±1.0 


2615 keV 


6.8±0.4 


6.6±0.6 


6.3±0.5 



3 a vs f3/j discrimination 

As described in Sec. [TJ the possibility to discriminate th< 
a interaction results to be the actual key point for a DBE 
bolometer. As described in details in [TH], in scintillat 
ing ZnMo0 4 bolometers the a vs /3/j discrimination can 
be obtained in two different ways: using the light signal 
and/or by using the PSD. We define discrimination power 
(DP) between the a and ,8/7 distributions the difference 
between the average values of the two distributions nor- 
malized to the square root of the quadratic sum of their 
widths: 



3.1 Light vs heat discrimination 

The light-to-heat energy ratio as a function of the hea 
energy is shown for the calibration spectrum in Fig. [3] /3/' 
and a decays give rise to very clear separate distributions 
In the upper band, ascribed to /J/7 events, the 2615 keV 7 
line is well visible. The lower band, populated by a decays 
shows the continuous background induced by the degradec 
a source. 

The evaluated LY of the ZnMo0 4 crystal is 1.54±0.0: 
keV/MeV. This value is surprisingly larger (+ 40%) witl 
respect to our previous measurements on 30 g (colorless 
samples [T^]. This, unexpected, larger LY results in ai 
increase of the DP using the scintillation light. Consider 
ing the events in the 2.5-^3.2 MeV region (see Fig. [3]) we 
estimate a DP of « 19. Moreover the scintillation yield 
evaluated on the internal a-line of 210 Po (see Sec. 0]) is 
0.257±0.002 keV/MeV that corresponds to a scintillation 
Quenching Factor of 0.167±0.002. This value is fully com- 
patible with the one (0.18±0.02) obtained on the 30 g 
sample. 

3.2 Pulse Shape Discrimination 

The PSD performed in this work is obtained with the same 
method described in [T^]. In Fig. 0] we plot the TVR vari- 
able as a function of the energy for the same data sample 
of Fig[3] As for the case of Figj3l (3/j and a events are dis- 
tinctly separated. The obtained DP with the PSD is « 14, 

2 Since we attribute to the heat peaks the nominal energy 
of the calibration 7's, the light-to-heat energy ratio also repre- 
sents the Light Yield of the crystal. 




3500 4000 



Energy [keV] 

Fig. 3. The light-to-heat energy ratio as a function of the heat 
energy obtained in the 80 h 228 Th calibration with ZnMo0 4 - 
Sum. The upper band (ascribed to /3/7 events) and lower band 
(populated by a decays) are clearly separated. The 2615 keV 
208 Tl 7- line is well visible in the /3/j band as well as a the con- 
tinuous background induced by the degraded a source. The 
events belonging to the energy region 2.5-^3.2 MeV (high- 
lighted in the plot) are used to evaluate the DP, that results 
w 19. 




ammm m 

500 1000 1500 2000 2500 3000 3500 4000 

Energy [keV] 

Fig. 4. TVR as a function of the energy, for the same events of 
Fig. [3] The upper band is populated by a particles (events in 
the 2.5-^3.2 MeV energy range are shown in red) while fi/^/'s 
contribute to the lower band (events in the 2.5-^3.2 MeV energy 
range are shown in black). 



as shown in Fig. [5) It has to be pointed out the "opposite" 
behaviour of the discrimination power with respect to the 
results obtained on the 30 g small samples [13] in which 
we obtained « 8 with the light signal and « 20 with the 
PSD. This is due to two distinct mechanisms. First the 
crystal tested in this work (despite its larger size and or- 
ange tint) emits « 40 % more light with respect to the 
small (colorless) sample previously tested. This implies, 
obviously, an improved DP of the scintillation light. Sec- 
ond the PSD is sensitive to the S/N ratio. In this work 
the FWHMf,ase resolution is 4 times worse with respect to 
the sample previously tested. 
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Fig. 5. TVR histogram of the events of Fig. [4] the a sample in 
red and the fi/j in black. The mean values and the standard 
deviations, as estimated from a Gaussian fit, are reported. The 
ft/j distribution was fit with an asimmetric Gaussian. The 
sigma of the right part of the gaussian (erf )is used to evaluate 
the DP. From Eq.[T]a discrimination power of « 14 is obtained. 

4 Internal contaminations 

The internal radioactive contaminations of this crystal 
were evaluated summing up background and different cal- 
ibration runs for a total collected statistics of 524 h. The 
corresponding a-spectrum is presented in Fig. El We found 
a contamination of 226 Ra (it shows a very clear a and 
"BiPo" decay pattern sequence). This contamination is 
evaluated as 27±6 /iBq/kg. As often happens, we also 
found a clear internal contamination of 210 Po, correspond- 
ing to an activity of 700±30 /iBq/kg. No other a lines 
appear in the spectrum. 

In order to evaluate the limits on other potential dan- 
gerous nuclei (in particular the 232 Th chain), we evaluated 
first the flat a continuum in an energy region in which no 
peaks are expected (3.6-^4 and 4.35-^4.7 MeV). Then we 
studied an interval of ±3ct centered around the Q value 
of each possible radioactive nucleus, being a the energy 
resolution of the 210 Po peak (4.5 keV). The expected fiat 
background contribution in each 27 keV energy window 
is 0.68 counts. Applying the Feldman- Cousin method [29] 
using the observed number of counts in each energy win- 
dow with the expected background, we were able to set 
90% CL limits on several nuclei, as reported in Tab. [3] 

5 Conclusions 

For the first time a large mass ZnMo04 crystal was tested 
as a scintillating bolometer for a possible next generation 
neutrinoless double beta decay experiment. The bolometer 
shows an excellent energy resolution. We demonstrated 
that, even on large mass detector, this compound is able 
to discriminate a particles interactions at -practically- any 
desirable level, using the light information as well as the 
pulse shape discrimination alone. Moreover this crystal 
shows an excellent radiopurity. 



> 




Energy [keV] 

Fig. 6. a-spectrum obtained in 524 h of measurement. The 
contribution of the 238 u/ 234 U a source is clearly evident below 
4 MeV. The internal a-lines arising from 226 Ra decay chain are 
highlighted. 

6 Acknowledgements 

This work was made in the frame of the LUCIFER ex- 
periment, funded by the European Research Council un- 
der the European Union's Seventh Framework Programme 
(FP7/2007-2013)/ERC grant agreement n. 247115. Thanks 
are due to the LNGS mechanical workshop and in particu- 
lar to E. Tatananni, A. Rotilio, A. Corsi, and B. Romualdi 
for continuous and constructive help in the overall set-up 
construction. Finally, we are especially grateful to Maur- 
izio Perego for his invaluable help. 



References 

1. A. Strumia and F. Vissani, |he~p-ph/0606054V 3 (2010). 

2. F.T. Avignone III, S.R. Elliott, J. Engel, Rev. Mod. Phys. 
80 , 481 (2008). 

3. W. Rodejohann, Int. J. Mod. Phys. E 20, 1833 (2011). 

4. A.S. Barabash, Phys. of Part, and Nucl. 42, 613 (2011). 

5. M. Auger et al, I arXiv: 12057 5608 vl (2012), submitted to 
PRL. 

6. A. Gando et al, Phys. Rev. C 85, 045504 (2012). 

7. S. Pirro, Eur. Phys. J. 27, 25 (2006). 

8. L. Cardani et al, JINST 7 P01020 (2012). 



Chain nuclide activity 



/iBq/kg 




232rj-,^ 


< 8 




228^ 


< 6 


238 jj 


238 Tj 


< 6 




234 u 


< 11 




230^ 


< 6 
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